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a^ . 

0^ , We investigate the resonant interaction of neutrinos in matter in the presence of mixing. 

At energies near the W resonance, oscillations are absent; the survival probability of 
electron antineutrinos is suppressed due to the presence of open inelastic channels. 



The existence of a resonant interaction of electron antineutrinos in matter and 
the possibility of flavor oscillations have been known for a long time; for a summary 
see [|l]. The laboratory energy at which the resonance occurs is about 6.4 PeV. 
^ I Recently, interest has turned towards an investigation of neutrino oscillations in 

^ I the PeV energy range, at energies near the W resonance in the interaction cross 

■ section @ . Learned and Pakvasa in ref . point out that neutrino oscillations in 

a hitherto unexplored range of neutrino masses and mixing angles may lead to the 
arrival and detection of r neutrinos in present and planned neutrino telescopes 
(DUMAND, NESTOR, KM3 etc. ). 

It is generally assumed that neutrinos emerging from a celestial source 



o 



p (typically, an AGN or a binary system) are produced as a result of the decay of 

p • hadrons. If so, in either hadron-hadron interactions or in photoproduction most 

D . of the hadrons produced are pions, with charged pions of either sign being pro- 

^ ' duced in approximately equal numbers and energy distributions. Consequently, 
the number of neutrinos and antineutrinos of any flavor is approximately equal. 

^ ■ To a reasonable approximation, the ratio of neutrinos (antineutrinos, resp.) of 

■ different flavors produced is: 

It has to be noted, however that at least some of the time, a considerable 
amount of matter is present between the primary source of neutrinos and the 
detectors in the experiments mentioned. 

• In all probability, neutrinos produced in an AGN or a binary system, have 
to pass through a region of high electron (and perhaps baryon) density. 

• In addition, the neutrinos usually pass through a certain amount of matter 
between the source and the detector, e.g. the Earth (in the case of the 
observation of upward going neutrinos) and/or sea water (in case of zenith 
angles close to the horizon). 



^Permanent address. E-mail: SKD@JHUP.PHA.JHU.EDU 
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In this note we begin to investigate the question of neutrino oscillations if 
neutrinos in the PeV energy range (around the energy of where the W resonance 
occurs) penetrate matter between the source and the detector. Our main purpose 
here is to explore the qualitative aspects of the phenomena. For this reason, we use 
the formalism developed by Wolfenstein and a hypothetical situation of having 
two flavors only. Furthermore, the medium is assumed to be a homogeneous one. 

A full treatment taking into account all three flavors and a more complete 
treatment of the presence of inelastic channels and inhomogeneities in the medium 
will be given elsewhere. 

The Wolfenstein equation can be viewed as a Schrodinger equation in a null 
plane formalism. (It has been known for a long time that the kinematics in the two 
dimensional transverse space is Galilean; for a description of the formalism with 
an emphasis on its group theoretical aspects, see ^, particularly the second of 
these references.) 

The behavior of the evolution operator of neutrinos, E, is governed by the 
equation: 

idtE = ^ ^E, 

2p 
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where p and P stand for the longitudinal and transverse components of the mo- 
mentum, respectively; and V are the mass squared and interaction operators. 
The variable t may be thought of as the time conjugate to the Hamiltonian, 

p 2 ^ M2 + 

ti = . 

The evolution operator as well as the mass squared, interaction and trans- 
verse momentum operators are to be thought of as / x / dimensional matrices, 
where / is the number of flavors. The operator P is proportional to the unit 
matrix. The important feature of the null plane formalism is that the longitudi- 
nal and transverse degrees of freedom are almost exactly separated. Thus, p may 
be diagonalized separately and quantities like V ^ governing the dynamics in the 
transverse plane may depend on its eigenvalue. (For the same reason, p is, to a 
good approximation, equal to the kinetic energy of the neutrinos at the source.) 

In the conventional treatment of neutrino oscillations (for a review, see [0]) one 
omits the square of the transverse momentum and represents the matter present 
in terms of a potential independent of the transverse degrees of freedom. (This 
is equivalent to the classic Lorentz formula, expressing the polarizability of a 
medium in terms of the forward scattering amplitude of light on the particles in 
the medium.) 

Since we are to treat a resonant interaction, such an approximation is not an 
adequate one. We can still neglect the transverse degrees of freedom; however, a 
more accurate expression has to be found for the "potential". The appropriate 
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procedure has been known for a long time, see [§] and the original references 
quoted there. 

Briefly, the recipe is that wherever a potential appears in the the lowest order 
treatment of the scattering, it has to be replaced by an effective potential. The 
latter is given by the level shift operator (the "-ft" matrix") of the scattering of 
the propagating particle in the medium. The relationship between the transition 
operator, T and the K matrix is: 

T = -K(l+iK)"^ (2) 

We refer the reader to Brueckner's lecture in ref. [§] for a derivation of the 
"recipe"just quoted. However, we can make it plausible by recalling that the K 
matrix obeys the Lippmann-Schwinger equation. In an operator form the latter 
reads: 

K = V + V^-K, (3) 
e - Ho 

where e is the energy (in the present context, e = p) and Hq is the free Hamilto- 
nian, see e.g. ref. ||^. The symbol of the principal value, V, simply means that 
those eigenstates of the free Hamiltonian for which the inverse oie — Ho is singular, 
are to be deleted from a complete set. One now realizes that the first term in the 
iterative solution of eq. (H) is indeed the potential; moreover, in the absence of 
open inelastic channels, the operator K is a Hermitean one, thus it qualifies as a 
"potential" in a Schrodinger equation. 

Armed with this knowledge, we can now construct a phenomenological effective 
potential describing the resonant interaction of electron antineutrinos in matter. 
We describe the scattering amplitude (the matrix element of T) by means of a 
Breit-Wigner formula; the forward amplitude is given by: 



M2 - s - iUTt 



(4) 



In the last equation, M stands for the mass of the W, Ff is its total width. Gp 
stands for the conventional Fermi coupling constant and s is the square of the 
center of mass energy. In terms of the mass of the electron and p it is given 
by s = 2mf,p. We note that in terms of these quantities, the elastic width of 
the W is given by the expression, Fg = M^Gp/6tc\/2. In what follows, we take 
the numerator of eq. (^ at s = M^; this is a good approximation for a narrow 
resonance. Due to the fact that we neglect all but the resonant partial wave in 
eq. (^ , the expression of the matrix element of the K matrix is found immediately, 

VIZ. 

MFe 

M2-s-^MFi' 

where Fj = F^ — Fe is the inelastic width of the W. If the inelastic width, Fj = 
Ft — Fe were zero, the expression in eq. (|) would be real as required. Thus the 
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effective potential in the medium entering eq. (ffl) is given by the expression, 



rip 



me 



(6) 



where is the electron density in the medium. 

In the present model, the ffavor space is two dimensional; we denote the ffavors 
by e and /i, respectively. Due to the fact that (in real life), both the electron and 
muon are much lighter than the W, the matrix elements of the effective potential, 
T are: 

e) = (e|J^| e) = F (7) 

and all other matrix elements vanish. 

Next, we give the matrix elements of the mass squared matrix, A^^: 

{e\M''\e) = -{fi\M^\fi) 

= cos 29 

{fi\M^\e) = {e\M^\e) 

= sin 26 (8) 

As usual, we subtracted the average mass squared; its presence merely multiplies 
the state vectors by an overall phase factor, hence, it has no influence on neutrino 
oscillations. 

We now have to solve eq. (|I]) with the effective Hamiltonian put together 
from eqs. and (^. Despite the fact that the Hamiltonian is not a Hermitean 
matrix in flavor space, the equation can be solved immediately. This is due to 
the fact that any non singular 2x2 matrix can be diagonalized by means of a 
similarity transformation. (If the matrix in question is not a Hermitean one, the 
diagonalizing matrices are elements of the group SL(2,C) instead of SU(2).) We 
merely quote the result of the calculation. 



(e I -E| e) = exp —it 



4p 



cos — z sm ( 



A^ cos2e + F/2' 



(u \E\ mu) = exp —it — 



,A2cos2e + F/2 
cos (p + zsm (p 



, , . f . F\ . A2sin2e + F 
(u Lfc e) = — zexp —tt — sm<p 

, , f ,A2sin2e 
(e \E\ u) = — zexp —it — sm0 

The quantities /i^ and entering these equations are defined as follows. 



(9) 



/i^ = (A^ + FA2(cos2e + sin2e) +FV4 

= ^ 
^ 2p 



1/2 
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The physical interpretation of these equations is not quite easy, due to the fact 
that the quantity F is complex. However, near the resonance itself (where the use 
of the Breit-Wigner formula is justified in the first place), considerable simpli- 
fications occur. In order to understand this, we notice that the quantity n^jnif, 
ocurring in the definition of the effective potential, defines a mass scale charac- 
terizing the medium. In the following Table we exhibit approximate values of the 
mean electron density, together with the characteristic mass scale, m^ = n^jme 
for some environments of interest. 



Electron densities and characteristic masses 
for some environments 



Environment 


ne 


cm ^ 




m,^[eV^] 


stellar interior (sun) 
Earth 
water 


1.6 X 10^4 
3 X 10^3 


2 X 10^] 

3 X lO'^ 
5 X 10^ 



We see that the characteristic mass scales are, in general, considerably larger 
than the expected values of in a typical mixing scheme occurring in the liter- 
ature, see, for instance, |]l| for a review. As a consequence, the expression of /x^ is 
considerably simplified near s = M^. In fact, we have: 

fi^ ^^ + A^ (sin 29 + cos 29) , (10) 

so that to leading order, /i^ ^ F/2. Using the Euler decomposition of the 
trigonometric functions, one realizes that in the terms proportional to 

exp {—itF/ip) expicj), 

F cancels out to a good approximation. By contrast, in terms proportional to 

exp {—itF/Ap) exp — i0, 

F is present in the exponent. The quantity F has an imaginary part due to the fact 
that there are inelastic channels open: W decays mostly hadronically. Therefore 
terms proportional to exp{—iFt/2p) are exponentially damped and all surviving 
terms in the expression of the transition matrix have the same phase. (By putting 
in numbers one concludes that the extinction distance is a few metres at energies 
near the resonance in any environment of interest; thus, at any reasonable value of 

the damped terms are negligible.) Consequently, we omit terms proportional 
to exp(— ImFt) in what follows. 

We now use the approximation introduced in eq. (p!OD to give the expressions of 
the transition probabilities. This is best done by introducing the dimensionless dis- 
tance in energy from the resonance, x and an appropriately scaled ratio of A^/m^: 
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The transition probabilities are then given by the expressions: 

(5 sin 20)^ 



\{e\E\e)[ 
mE\e)[ 



1 + (x + 25 (sin 20 + cos 26) (1 + x'^)Y 
l + (a; + (5sin2e(l+a;2))^ 

1 + (x + 2(5 (sin 26 + cos 26) (1 + x'^)f 
l + (a; + 5sin26(l+x2))^ 

1 + (x + 25 (sin 26 + cos 26) (1 + x'^)f 



(5sin26)^(l + x2 



\{e\E\ii)\' = ^ '- 2 (12 

1 + (x + 25 (sin 26 + cos 26) (1 + a;2))^ 

By looking at the last set of equations, one realizes that the survival probability 
of electron antineutrinos is very low; this is due to the fact that the W resonance is 
a highly inelastic one: once the electron antineutrino interacts, the W decays into 
the original channel only rarely. By contrast, the muon antineutrinos are hardly 
depleted and they hardly ever produce electron antineutrinos: these processes 
proceed through mixing only. Remarkably, the t-dependence of all transition 
probabilities disappears beyond the extinction distance: an inelastic resonance 
produces uniform transition probabilities. 

Let us summarize. 

• The presence of an inelastic resonance stongly suppresses oscillation effects. 

• The original 1:2:0 ratio of the antineutrinos of the three flavors is strongly 
distorted. Electron antineutrinos are depleted, while muon antineutrinos 
are hardly affected. (In fact, their number is slightly increased due to the 
decay W — >• iip^.) 

• One can reasonably conjecture that in a calculation taking three generations 
into account, Vr production will be comparable to the rate (with the 
appropriate mixing angles and mass differences), i.e. a considerable number 
of T leptons and the corresponding antineutrinos will be produced. In par- 
ticular, the appearence of Pt--s in a neutrino telescope does not necessarily 
imply oscillation lengths comparable to the distance to some AGN. 

• The ratios of neutrinos (as opposed to antineutrinos) should be unaffected 
unless there exist neutrino-antineutrino oscillations. 
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